Singlet and Triplet Superfluid Competition in a Mixture of Two-Component Fermi 

and One-Component Dipolar Bose gases 
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We consider a mixture of two-component Fermi and (one-component) dipolar Bose gases in which 
both dipolar interaction and s-wave scattering between fermions of opposite spins are tunable. We 
show that in the long wavelength limit, the anisotropy in the Fermi-Fermi interaction induced by 
phonons of the dipolar condensate can strongly enhance the scattering in the triplet channel. We 
investigate in detail the conditions for achieving optimal critical temperature at which the triplet 
superfluid begins to compete with the singlet superfluid. 



The ability to easily mix cold atoms of different species 
to form new quantum systems brings another exciting di- 
mension to the study of ultracold atomic gases. A single- 
component Fermi gas only supports Cooper pairing with 
odd parities, such as p-wave pairing, which are typically 
strongly suppressed in accordance with Wigner's thresh- 
old law p]]. Mixing bosons induces an attractive interac- 
tion between fermions |2j , which, as Efremov and Viverit 
Q pointed out, raises the prospect of achieving p-wave 
superfluidity in a Fermi-Bose (FB) mixture. Recently, 
Dutta and Lewenstein [4| generalized the idea to a (2D) 
mixture involving dipolar bosons with the goal of real- 
izing a superfluid with p x + ip y symmetry whose excita- 
tions are non-Abelian anyons that are the building blocks 
for topological quantum computation [5| , and Nishida [6] 
sought the same goal by mixing fermion gases of different 
species in different dimensions. 

In this Letter, we consider a (3D) homogeneous 
mixture (with an effective volume V) between a two- 
component Fermi gas and a dipolar Bose gas, made up 
of two equally populated (balanced) hyperfine spin states 
( | f) and |4_)) of a non-dipolar fermionic atom of mass m p , 
and a ground state of a bosonic molecule (or atom) of 
mass rriB with an induced dipole aligned along the ex- 
ternal electric field direction z. The two pseudo spins 
provide fermions with the opportunity to pair not only 
via triplet (with odd parities) but also via singlet (with 
even parities) channels of interaction. This opens up the 
possibility of using this two-component FB model to em- 
ulate and explore pairing physics analogous to that in 
superfluid 3 He [7J, which is known to be greatly enriched 
by the existence of an internal spin degree of freedom. 

In the context of ultracold atomic physics, there has 
been a recent upsurge of activity in pursuing similar goals 
but with (3D) two-component dipolar Fermi gases [f§, 
motivated largely by recent rapid experimental advance- 
ment in achieving ultracold dipolar gases both in 40 K 
- 87 Rb polar molecules [l(|, and in Cr 11] and spin-1 
Rb atoms 12]. Such studies [j, Q represent a gener- 
alization of earlier work 13j, aimed at exploiting a 
rf r 2_ 3z 2 -type of anisotropy in dipole-dipole interactions 
for enhancing triplet pairing in single-component dipolar 
Fermi systems. 



The induced Fermi-Fermi interaction mediated by a 
dipolar condensate is also anisotropic in nature and thus 
opens up a new avenue for studying superfluids with 
unusual pairings. The progress in this area has so far 
been limited, to the best of our knowledge, to a single- 
component model in a 2D geometric setting Q. In con- 
trast, the present work expands such studies to a 3D 
two-component mixture, where both the dipolar interac- 
tion between bosons and the s-wave scattering between 
fermions of opposite spins are independently tunable, and 
seeks to use it as a model to explore the physics that are 
currently being hotly pursued in two-component dipolar 
Fermi gas systems |9|. In this Letter, we study in detail 
the anisotropic nature of the 3D induced interaction, and 
particularly how one should prepare a two-component FB 
mixture in order to maximize the opportunity this in- 
duced interaction affords for raising critical temperatures 
at which phases of different parities begin to compete. 

To begin with, we model our system with a (grand 
canonical) Hamiltonian H — Hb + Hbf + Hp: 

Hb = J2 = e ^ B - Vb) blk + (2V)- 1 x 
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where 6 k (o-o-.k) is the field operator for annihilating a 
boson (a fermion of spin a) of kinetic energy e k ,_B = 
H 2 k 2 /2mB (ek,F = fi- 2 fc 2 /2mp) and chemical potential 
Mb (a*f)- The low temperature physics of the mixture 
under consideration arises from the interplay between 
short- and long-range two-body interactions. The for- 
mer is dominated by s-wave scattering characterized with 
strengths: Ubb — 4:TtH &b b/wibi Ufp = Aith 2 a,FF fmp, 
Ubf = 4:irh 2 aBF l"mBp\^ 2mB"m,p / '(tub + ^f)], and 
U<j<t' (q) = U ff^ct',— cr? where clbb, &bf and app are 
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the related scattering lengths, and the Kronecker-<5 func- 
tion in U aa > (q) limits the fermionic s-wave interac- 
tions to fermions of opposite spins. The latter is the 
dipole-dipole interaction (restricted to bosons) given by 
Udd (q) = 8nd 2 P2 (cos# q ) /3 in momentum space, with 
d 2 the dipolar interaction strength, P 2 (x) = (3x 2 — l) /2 
the second-order Legendre polynomial, and # q (0 q ) the 
polar (azimuthal) angle of vector q. 

In the low temperature regime considered in the 
present Letter, bosons are virtually all condensed to 
the zero-momentum mode and a straightforward appli- 
cation of the Bogoliubov approximation, in which &k=o is 
treated as a c- number 6k=0i yields a well-known picture of 
the bosonic system described by Eq. ((T|): it consists of a 
collection of phonon modes with the Bogoliubov disper- 
sion relation E k = VB^k\Jl + (£sfc) 2 + 2sdd,P2 (cos 9 k ) 
(lij and a homogeneous dipolar condensate with den- 
sity riB — |&k=o| 2 (^b = n B Ub b ) j which is stable 
against collapse provided £dd < 1, above which phonons 
with k — > acquire imaginary frequencies. Here, 
£dd = And 2 /(3Ubb) measures the strength of the 
dipolar interaction relative to the s-wave interaction, 
Vb = \/tib Ub b/ttib is the phonon speed, and £s = 
h/ ^/AmsnsUBB is the healing length. Integrating away 
the phonon degrees of freedom [2j leads to an effective 
Fermi system described by the same Hamiltonian as Eq. 
([3]), except that Uaa> (k) = Uff^'-g + U in d (k), where 



Uind (k) 



U B f/Ubb 



1 + (6j*0 +2£ dd P 2 (cos0 k ) 



(4) 



is the phonon-induced Fermi-Fermi interaction in the 
static limit @, Q . As can be seen, the induced interaction 
depends on the dipole orientation differently than the di- 
rect dipole-dipole interaction and therefore provides an 
alternative model for the exploration of spin singlet and 
triplet paring. 

Typical BCS mean-field theory proceeds with the in- 
troduction of the matrix representation for the BCS or- 
der parameter in the uncoupled spin space: A a i a (k) = 
J2w Ucrcr> (k — k') (<Z-k',CT'Sk',cr)- The notion of spin sin- 
glet and triplet pairing emerges when one moves from un- 
coupled to coupled spin space spanned by a spin singlet 
\S = 0, M = 0) and triplet \S = 1, M = -1, 0, +1), which 
are antisymmetric and symmetric with respect to the spin 
exchange, respectively, where M is the z projection of the 
total spin S. As one may easily verify, the gap parame- 
ter, A s (k) = (A n (k) - A it (k)) /2 associated with the 
singlet (S — 0) state is an even function of k, while the 
three gap parameters, A l ' x (k) = (A^ (k) — A-^ (k)) /2, 
A*.»(k) = (A u (k) + A n (k))/2i, and A^ (k) = 
(A-fj, (k) + Aj/i- (k)) /2, associated with the triplet (S — 
1) states are odd functions of k, in accordance with 
Fermi statistics, where use of A a p (k) = A s (k) i (<J y ) a p + 
J2u=x. y .z A*'" (k) i (<T u a y ) a p, a convention in the study 



of superfiuid 3 He 0], has been made, with a u being the 
usual Pauli matrices. 

To highlight the dominant physics, we ignore the Fermi 
surface deformation due to the anisotropy of the Fermi- 
Fermi interaction 17J and consider two-body scattering 
up to the level of the Born approximation 1J, 14 1 , both 
of which hold in the weakly interacting regime where 
tifU^p /Ubb <C c_f = Fermi energy. At critical tem- 
peratures where the gaps are small, one can ignore the 
nonlinear coupling between parings of different parities 
so that the critical temperatures can be estimated by a 
set of linearly coupled gap equations 0, : 



A (k) = -V~ l U (k, k') K {k') A (k) , 



(5) 



with the understanding that A (k) = A s (k) and 
[/(k, k') = J7 s (k, k') for singlet pairing and A (k) = 
A*' u (k) and U (k, k') = IP (k, k') for triplet pairing, 
where U s (k, k') = U FF + [U ind (k - k') + U md (k + k')]/2 
and U* (k, k') = [U md (k - k') - U md (k + k')] /2 are the 
singlet and triplet potentials that are even and odd 
functions of both k and k', respectively, and K (k) = 
tanh (/3&/2) / (2&) - l/(2e fc ). In arriving at Eq. ©, 
we have applied the standard procedure to renormalizc 
the contact interaction [2j and a similar procedure (but 
expressed in terms of vertex functions [14J) to renormal- 
ize the dipolar interaction. 

Making decompositions: A (k) = J2i A; (k) 5^°(k) and 
L/(k,k') = ^T.i,v, m Ui m>Vm {k,k')Y{ m {^)Y Vm {ii'), in 
which the azimuthal symmetry of the interaction ((4]) has 
been explicitly incorporated, we cast Eq.([5]) into the form 
containing only the radial coordinate: 

A * ( fc ) = - E / ^ K ( fc/ ) u w° ( fc > Ai ' 

(6) 

where Yj m (k) are spherical harmonic functions. In the 
low temperature limit ksT / ep <C 1, k 2 K (k) is small vir- 
tually everywhere except around the Fermi momentum 

1/3 

kp [= (37r 2 n^) ] where it is sharply peaked compared 
to other momentum distributions, and the critical tem- 
perature can be estimated, to a good approximation, by 
the equation for the gap parameter, A; = A; (kp), at the 
Fermi surface 

nksT 



Ai = N (cf) In ■ 
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Uu.Ai 



where 7 = 0.577 is Euler's constant, and Ui y = 
Uioj'o (k F , kp) is an element of the interaction matrix 



U: Uf >v = U F F5ifiSv fi + 
0.5[1 - (-lfpffi, with 



U u , = 2tt 



0.5[1 + {-l) l ']Ujf and UK = 



U ln d(k F k - k F k')d((j) k - 0k') x 
Yio (cos6» k ) Y/'o (cos 9 k >)d (cos 6 k )d (cos 6v) . (8) 
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FIG. 1. The scaled eigenvalue uj 3 for the singlet state (a) and 
uj t for the triplet state (b) as functions of 5 for different Edd- 
In both figures, the dotted curves are for the nonpolar case 
(edd = 0) and the solid curves from bottom to top are for (a) 
e d d = 0.6, 0.8 and 0.9, and (b) E dd =0.1, 0.5 and 0.9. 



In contrast to the dipole-dipole interaction, which couples 
A; only to A/±2, the dipole induced interaction in Eq. 
(|4]) can, in principle, couple A; to any A(; +2 n)>o with n 
being an integer according to Eq. (|8|). As expected, the 
singlet and triplet pairings are decoupled, containing all 
the even and odd partial wave components, respectively. 

To solve Eq. (J7|), we move into a primed space in 
which the interaction matrix U' — MUM* is diagonal- 
ized via a unitary transformation M so that U' , — 
—uj n {Ub F /Ubb) 5 nyn ' , where ui n is the eigenvalue scaled 
to — U bf /Ubb- In this primed space, Eqs. (J7]) are decou- 
pled, leading, immediately, to the critical temperature 



Tn 



it F e 



j-2 



■ exp 



uj n (6) A 



(9) 



for the nth channel, in which the order parameter 
A (k) oc J^i M* ;T;o(k) becomes a superposition of differ- 
ent partial waves with different angular momenta, where 



5 = £,BkF = k F / (4y^msOBs) , 

A = 2^)^ = 1™^ 
Urr tt m„„ as 



— ^k F . 

7r m BF a B B 



(10) 
(11) 



As the temperature is lowered, the most favorable su- 
pcrfiuid phases to be realized correspond to those chan- 
nels with the strongest attractive interactions (the high- 
est positive u> n ). 

Figure [T] shows how the strongest (attractive) inter- 
actions in the singlet and triplet channels, lo s and uit, 
change with 5 for different Edd when Uff = 0. The 
dotted curves, ui s = (26) ln[l + (26) ] and Ut — 
2 (26y 2 |ln[l + (26) 2 }[(26y 2 + - lj, represent the 
corresponding interactions in mixtures with nonpolar 
molecules (edd — 0) 0, HH, where, as 6 reduces, the 
triplet interaction begins to decrease to zero after reach- 
ing its peak around <5 w 1 as opposed to the singlet inter- 
action, which increases monotonically. The introduction 



of a dipolar condensate to a fermion gas adds to the de- 
nominator of Eq. ((4]) an anisotropic term 2eddPi (cos #k), 
which plays an increasingly important role compared to 
the isotropic contribution 1 + 6 2 in the small 5 region, 
around which significant changes are observed to take 
place in Fig. [1] The most striking development hap- 
pens, however, in the triplet interaction [Fig. [TJb)] which 
not only asymptotes to a finite value in the limit 5 — > 0, 
in clear defiance of Wigner's threshold law, but also in- 
creases with Edd for a given 6 far more dramatically than 
the singlet interaction [Fig. QJa)]. This provides concrete 
evidence that the use of a dipolar BEC can indeed signifi- 
cantly enhance scattering in the triplet channel compared 
to the singlet channel, which has been the key motiva- 
tion behind the present proposal for achieving the triplet 
superfluid. 

To prepare a system with small 6, we must employ a 
Bose component with a relatively small (high) healing 
length (density). This often means that the system sep- 
arates into a mixed phase with densities (riFi,nBi) and a 
pure Fermi phase with densities (n F 2, nB2 — 0), the only 
phase separation scenario that involves a mixed phase 
lil 20j |. (A complete separation between fermions and 
bosons requires much higher densities than considered in 
the present work.) The mixed phase must share the 
same chemical and thermodynamical potentials with the 
pure phase. This consideration leads to 



Ubfubi 

-U B Bn B i/^ - U BF n Bi n Fi - 2An 5 Ji/5 



A 2/3 i 2/3 



-2At44 3 /5 



from which one finds 



n B i 



An 2 ^(y'-1)/U 



BF, 



(12) 



where A — (37r 2 ) 2 ^ 3 h 2 / (2m F ) and y — (np2/npi) 1 ^ 3 is 
the solution to the cubic equation 

- 15 (y + l) 2 /A + 8y 3 + 16y 2 + 24y +12 = 0, (13) 

(see Ref. [lj| for details) . All previously derived formulas 
concerning the critical temperature are directly applica- 
ble to the mixed phase upon substitution of ub with ubi 
and np with np\. 

The optimal triplet superfluid temperature T t is always 
found to occur in the mixed state. An example in which 
mp = 6u, tub = 127u, and asB = 250ao (with u the 
atomic mass and ao the Bohr radius) is given in Fig. [2] 
(a), which illustrates how T t and the required asF change 
with npi for different Sdd- In arriving at Fig. [5] (a), we 
used Eq. ([9]) to construct, for a given set of rip\ and 
Edd, temperature T t as a function of A by solving for 
the required clbf, nsi, and 6 simultaneously from Eqs. 
(fTO]) - (fT3|) . The optimal T t corresponds to the peak 
temperature at some A = A pca k- In Fig. H](a), we find 
that Apoak ~ 0.6 and <5 pca k ~ 0.25 with small variation for 
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n Fi (cm" 3 ) U FF /(.Ujj F /U BB ) 

FIG. 2. (a) The optimal T t and the required hbf as functions 
of uf 1 for different edd under conditions that tuf — 6u, ms = 
127m and ass = 250ao. (b) illustrates how T s (curves) can be 
made degenerate with Tt (horizontal lines) by changing Uff 
while fixing all the other parameters to those represented by 
the black dots in (a) where risi have reached 5 x 10 14 cm" 3 . 
In both figures, edd = 0.9 for solid curves, edd = 0.85 for 
dotted curves, and edd = 0.8 for dashed curves. 



triplet superfiuids in a two-component FB mixture with a 
dipolar condensate. We have found that T s can be made 
degenerate to T t at a temperature 10 7 orders of mag- 
nitude higher than 10~ 6 nK (not shown), the optimal 
temperature achievable under a similar set of fixed pa- 
rameters for a two-component FB mixture with nondipo- 
lar bosons. Just as mixing nonlinear waves has been an 
important means for creating coherent sources of laser 
light, mixing cold atoms is expected to play an increas- 
ingly more important role in creating new quantum gases 
(or liquids) in the coming years as the field of ultracold 
atomic physics continues to mature. The present study 
reinforces the notion that mixing fermions with dipolar 
bosons adds another exciting dimension in the pursuit of 
quantum systems with new and novel properties. 

H. Y. L. acknowledges the support from the US Na- 
tional Science Foundation and the US Army Research 
Office. 



different values of npi and edd- As can be seen from the 
solid curves (edd — 0.9), a temperature about 15 nK can 
be achieved in a mixed phase (marked with a black dot) 
with densities (n B i = 5 x 10 14 cm" 3 , n F1 = 2.17 x 10 12 
cm~ 3 ) and clbf — 304ao- 

As to the singlet supcrfluid temperature T s , it depends 
on the contact interaction Up p [or (kpcipp )~ ], which in 
our model is made magnetically tunable via Feshbach 
resonance. Thus, in principle, the interaction in the sin- 
glet channel can be made degenerate to that in the triplet 
channel by tuning both the dipolar interaction (edd) with 
an electric field and the s-wave scattering length (ciff) 
with a magnetic field. This opens up the possibility of 
studying phase competition between singlet and triplet 
superfiuids, a recurring theme in current studies con- 
cerning two-component dipolar Fermi gases. Figure [5] 
(b) illustrates how T s (curves) can be made to cross T t 
(horizontal lines) by changing Upp for different Edd- I n 
contrast to the pure two-component dipolar Fermi gas 
model, where, due to the average of the dipolar interac- 
tion over all the directions being zero, the singlet interac- 
tion is always less attractive than the triplet interaction 
in the absence of Uff, and Uff must be tuned to the 
negative side of the Feshbach resonance in order to make 
T s comparable to Tt [9J], the singlet interaction in our 
model is more attractive than the triplet interaction in 
the absence of Uf f [Fig- Q] , and consequently only when 
Uff is tuned to the positive side of the Feshbach reso- 
nance, can T s be brought down to the level of T t [Fig. [5] 
(b)]. An estimate based on Uff/ (U bf /Ubb) ~ 1-08 at 
the crossing of the two solid lines (edd = 0.9) in Fig. [2] 
(b) indicates that T s = Tt when (kpap) 1 = 2.03. 

In summary, we have investigated the optimal condi- 
tions for achieving the coexistence between singlet and 
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